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Abstract 

In this paper, we examine the effectiveness of absorbing layers as non-reflecting com- 
putational boundaries for the Euler equations. The absorbing-layer equations are simply 
obtained by splitting the governing equations in the coordinate directions and introducing 
absorption coefficients in each split equation. This methodology is similar to that used by 
Berenger for the numerical solutions of Maxwell’s equations. Specifically, we apply this 
methodology to three physical problems shock-vortex interactions, a plane free shear flow 
and an axisymmetric jet with emphasis on acoustic wave propagation. Our numerical 
results indicate that the use of absorbing layers effectively minimizes numerical reflection in 
all three problems considered. 
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tract No. NAS 1-19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center. Hampton, VA 23681-0001 
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1. Introduction 

The proper treatment of computational boundaries is crucial for any numerical solution to 
a set of partial differential equations which governs fluid motion or wave propagation in a 
medium. Various techniques have been developed to minimize the reflection of out-going 
waves. A review can be found in Givoli (1991). Numerical boundary conditions based on 
the characteristics of the relevant linearized equations and their asymptotic solutions in the 
far field have been widely used. However, such boundary conditions are not satisfactory if 
the outflow is nonlinear or involves multi-directional waves. As a possible remedy, a buffer 
zone abutting the computational boundary, in which the governing equations are modified, 
and whose role is to absorb the incident waves, has been proposed. In this buffer zone, the 
modifications have the effect of either removing or damping reflected waves oriented back 
towards the computational domain. Naturally, the buffer zone solutions themselves need not 
necessarily be physical, and they serve only to prevent contamination of the solution in the 
physical domain of interest by the reflections from the computational boundaries. Various 
types of buffer zone techniques have been used in flow simulations. For example, Oolonius 
et al.(1993) used buffer zones in which the solutions were filtered. In a different approach, 
Ta’asan and Nark (1995) modified the governing equations in the buffer zone to change the 
orientation of the characteristics, and make the flow supersonic at the exit plane. Recently. 
Berenger( 1994) proposed a very effective Perfectly Matched Layer technique for Maxwell’s 
equations. In this approach, the equations governing the so-called matched layer are split 
into subcomponents with damping terms which absorb the incident waves almost perfectly. 
Following Berenger, Hu (1996), developed an analogous technique for the linearized Euler 
equations, and provided analytical results for the case of uniform flow. 

In this paper, we follow the operator splitting principle of Berenger (1994) and Hu (1996) 
for the equations governing what we call the absorbing layers and examine their effectiveness 
in the case of shock-vorticity wave interactions, a plane free-shear layer and an axisymmetric 
jet. The emphasis is on the effectiveness of the the computational boundaries in handling 
wave propagation including sound waves. It is shown that the absorbing layer technique is 
very effective for all three physical problems. The next section describes briefly the numerical 
models used in this study, followed by the section on results and conclusion. 


2. Numerical Models 

2.1 Shock Wave Interactions 

To verify the applicability of the absorbing boundary condition technique to shock-turbulence 
and shock-vortex interaction problems, we choose the numerical model of Erlebacher, Hus- 
saini and Shu (1997). This model solves the fully nonlinear compressible Euler equations 
along with a time evolution equation for the shock motion for the purpose of fitting the 
shock. The outflow boundary conditions which minimizes wave reflection back into the do- 
main of computation are of crucial importance for such problems as they involve long-time 
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integrations. The present case focuses on the interaction of a single vorticity wave with 
a shock wave, and the results of course carry over simply to a randomly distributed wave 
system. The two dimensional Euler equations are written as 


dp dp dp du dv , 

dt + U lh + '% = ~ P ^Ih + dy‘ 


du du du 1 dp 

dt dx dy p dr 


dv dv dv dp 

m +u »i + , % = ~ihi 


dp dp dp , du dv . 

Ii + u £ + % = - v{ aI- + Ty ) 

The computational domain has the shock as a boundary on the left and an outflow 
boundary on the right, and is periodic in the other direction. Fourth order Runge-Kutta 
scheme is used for time integration, and the spatial derivatives are discretized by a compact 
6 th order scheme. 

In the absorbing layer at the right boundary, the Euler equations are split into a lo- 
cally one-dimensional set with artificial damping terms. Consider the pressure equation, for 
example, in computational space: 

Op Op Op OiV] 0w 2 
_ = _ ai _^ ! __ TP( _ + _ ) 

After operator splitting and addition of damping terms, the pressure equation becomes 


dp i 

dp 
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and iv 2 are 
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components 


«i and a 2 are cont ravariant velocity components (which include the effect of grid motion) 
in computational space, and p = p\ + p 2 * Locally one-dimensional equations for the other 
variables are constructed in a similar manner. The damping factor ax is zero in a layer 
parallel to the A" direction; similarly ay is zero in a layer parallel to the V direction (see 
Figure 1). However, in the corner region both these damping factors are positive. 
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2.2 Free Shear Layer 


In order to evaluate the performance of the absorbing-layer technique in the case of inviscid 
instability waves, vve solve the linearized Euler equations in a Cartesian {x,y) coordinate 
system. We study the evolution of a Kelvin-Helmholtz instability wave as it propagates 
downstream and impinges on the absorbing layers. In this case, the ^-momentum equation 
reads 

du - du dU 1 Op 

777 + £' tt — h “T“? ? + ~tt~ = 0 
at ox dy p ox 


where 


U — ~ 2 [(C r i + ( 2 ) + (fq — U 2 ) tanh(y)]. Absorbing layers are used at the upper, lower and 
right boundaries. Again, the afore-mentioned operator splitting in the absorbing layer leads 
to two ^-momentum equations: 


dii\ -du 1 Op 

”77T “f &x u \ + t 1 77 h — TT~ — 0 

at dx p dx 


du 2 dU 

~m + ‘ T »" 2 + V’ ~ 


where a — u\ + u 2 * All other equations are treated similarly. These equations are solved 
by a low-dissipation and low-dispersion Runge-Kutta scheme which is formally fourth order 
accurate (Hu, Hussaini and Manthey, 1996). 


For the nonlinear case one uses again an approximate time independent mean flow to 
split the Euler equations in the absorbing layer. Thus the stream-wise velocity for two 
dimensional flow's is decomposed into three components: 


n = (I + U \ + u 2 


where U is the mean velocity as in the linear case. Then the .r-momentum equation is written 
as 

du 1 

— + uu x + vu y + -p x = 0. 

dt p 

This equation is then split into two equations as 

du\ 1 1 

- 7 — + VxUi - — p x - uu x + -P x + UU X 
at p p 
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All other equations in the absorbing layer are similarly derived. 
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2.3 Axisymmetric Jet 

The compressible axisymmetric Euler equations for the jet in the weak conservation form 
are : Q t + F z + G r = S', where, in the linearized case, 
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In the above equations, p, p, m z , m r , E denote the fluctuating components of pressure, 
density, axial and radial momentum, total energy and H is the mean enthalpy. These equa- 
tions have been linearized around the mean velocity (U r , U z ) represented by an error function 
that fits experimental measurements. The interior equations are simply split into 

Qt + F z = 0, Q t + G r = S' 

and they are modified in the absorbing layer as 

Ql +F:= ~<T : Q' + S 1 , Qt + Gr = -°rQi l + S ' 2 

where Q = Q 1 + Q 2 and S = S’* + S’ 2 . (We used S 1 = 0 in this study). 

In the nonlinear case, the vectors O. S, /' . and G are defined as follows. 
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We use the fourth-order MacCormack method which has been successfully used in earlier 
studies by Hayder et al. (1996) to solve the linearized Euler equations, and by Hayder et al. 
(1993) and Mankbadi et al. (1994) to solve the Navier Stokes equations. The equations are 
linearized before splitting to obtain the equations for the absorbing layer. Thus, we get 

(Q' - Ql)t + (F- F 0 )~ = -a z (Q l - Ql ) 

( Q 2 — Qo)t + ((• — Go)r = —&r(Q 2 — Qo) + ( S — So) 
where the subscript 0 denotes mean quantities. 

3. Results 

In the case of the shock-vorticity wave interaction, we consider specifically the following 
simple wave 

u — U\ = tU\2y/2-^cos(k r x + k y y - U t t) 
v = —tl r , 2 \pl -j- cos ( k T x + k y y — U\t) 

p = p = T = 1 

as the upstream condition ahead of the shock. Ui is the upstream mean velocity normal to 
the undisturbed shock, k y = k sinO , k x = k cosB (where k is upstream wavenumber), arid 
e = 0.001 measures the intensity of the w r ave. Our standard interior domain is 7.4 units long 
with 185 uniformly spaced grid points. We used 16 points on the coordinate axis parallel 
to the shock. An absorbing layer abuts the right outflow boundary. We introduce damping 
gradually in order to minimize any reflections due to the discretization in the absorbing layer. 
Unless otherwise mentioned, we use 0 = 30°, I\ = 2, and 25 grid points (= 1 unit in length) 
in the buffer layer for our computations. A snapshot of pressure in the interior domain at 
t — 20 is presented in Figure 2, which shows how well the out-going waves are absorbed with 
little reflection. To measure the contamination due to reflection, the solutions are compared 
with a reference solution obtained by computing the flow in a much larger domain with the 
same spatial and temporal resolution. We follow this methodology for all problems in this 
study. Figure 3 compares axial variation in pressure for two different size buffers against 
the large domain solution at t = 20. Because of modifications to the governing equations, 
the solution in the buffer layer is irrelevant. The solution in the interior domain for a buffer 
with 25 points is visually indistinguishable from the larger domain solution. In Figure 4, 
we show the rms error (E) in pressure at the ordinate four grid points upstream of the 
interface between the computational domain and the absorbing layer as a function of the 
layer thickness measured in the number of equidistant points. The error E is defined as 

F _ ioo /E; v =,(/> r -P ) 2 

iKnJV N 
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where p r is the pressure from the reference solution, \p r max \ is its maximum amplitude and N 
is the number of grid points in the y-direction ( N=16 in the current context). E measures 
the numerical error in the solution, which includes both direct and induced errors due to the 
interaction of residual reflections from the outflow boundary with the flow and the shock. 
As expected, E decreases as the layer width is increased. In Figures 5 and 6, we show the 
dependence of numerical errors on the angle of incidence (9) and the wave number (k). The 
buffer layer is more effective at lower incidence angle and wavenumbers, although we notice 
some cross-overs in our numerical experiments. At later times, a fraction of the reflections 
from the outflow' boundary propagates upstream. These waves can then reflect back and 
forth, and cause what we call induced errors. These sometimes constitute a significant 
portion of the errors shown in Figures 4-6 at later times. 

The results for the free-shear layer are obtained for upper and lower stream mean veloc- 
ities, normalized by the speed of sound, equal to l i = 0.6 and ( 'i — 0.2 respectively. The 
eigenfunctions of the Kelvin-Helmholtz instability wave given by the linear stability theory 
are forced at the inflow, with a maximum amplitude t equal to 0.01. We solve the linearized 
Euler equations and the solution agrees with the linear theory very well in eigenfunction 
and growth rate comparisons. In Figure 7, snapshots of axial velocity (Fig 7a.) and pressure 
(Fig 7b) are shown, and in Figure 8 we present the amount of reflection as a function of 
the layer thickness. We observe that for 10 points in the absorbing layer, the amount of 
reflection (measured four grid points away from the buffer layer boundary) is less than .03% 
of the amplitude of the reference pressure fluctuation from the large domain solution. We 
also solve nonlinear Euler equations where the nonlinearity in the flow is significant. The 
inflow excitation amplitude (f) is kept at 0.01, but the interior domain is three times longer. 
All other flow r parameters are the same as in the linear case. The error in pressure four grid 
points away from the buffer layer in shown in Figure 9. We needed a larger buffer layer for 
the nonlinear flow’ simulations. At time equal to 3000, errors w'ith 30 and 50 points in the 
buffer layer were 3.5% and 4%> respectively. Intuitively one expects that a buffer layer to be 
more effective if nonlinear effects are smaller. This may be the principal reason for larger 
errors in Figure 9 compared to Figure 8. The effect of nonlinearity is also shown in Figure 

10, where we compare errors for simulations with two different levels of excitation e w r ith a 
buffer of 50 grid points. 

Finally, for the case of the excited axisymmetric jet, we assume the mean Mach number 
to be 0.6. At the inflow, we extrapolated one characteristic variable corresponding to the 
outgoing acoustic w’ave from the interior and computed the other three characteristic vari- 
ables at time t using [ p,u, v,p ] = tRe{qe luJt ), where q = [p, u,v,p\ is the eigenfunction given 
by the linear stability theory, t = 10 -4 , u> = 1.05. A snapshot of pressure is shown in Figure 

11. The rms pressure error (E) in the immediate neighborhood (four points away from the 
buffer layer) of the layer interface is plotted in Figure 12 for time equal to upto 50. This 
error becomes quasi-periodic and the maximum error for 25 grid points in the absorbing 
layer is about 0.015%. Our results for the nonlinear Euler equations are shown in Figure 13. 
The domain size is 10 units long for both linearized and nonlinear Euler simulation of the 


6 



excited jet. The physical parameters are the same for both the linearized and the nonlinear 
Euler equations for the jet calculations. 


4. Conclusions 

In conclusion, we find the performance of the absorbing-layer technique in the cases of 
three physical problems (using three different numerical algorithms) is quite satisfactory. 
This technique offers a viable alternative to the traditional boundary treatments based on 
the linearized characteristics or asymptotic solutions in the far field, and also other types of 
buffer layers. It also promises to be accurate and inexpensive for aeroacoustic computations. 
Further studies are warranted to put this methodology on a firm footing. 


References 

Berenger, .J-P, “A Perfectly Matched Layer for the Absorption of Electro-magnetic Waves”, 
Journal of Computational Physics , 114: 185-200, 1994. 

Colonius, T., Lele, S. K. and Moin, P., “Boundary Conditions for Direct Computation of 
Aerodynamic Sound Generation”, AlAA Journal , 31: 1574-1582, 1993. 

Erlebacher, G., Hussaini, M. Y. and Shu, C., “Interaction of a Shock with a Longitudinal 
Vortex”, ICASE report 96-31, to appear in Journal of Fluid Mechanics , 1997. 

Givoli, D., “Non-reflecting Boundary Conditions”, Journal of Computational Physics. 94: 
1-29, 1991. 

Hayder, M. E. Turkel, E. and Mankbadi, R. R., “Numerical Simulations of a High Mach 
Number Jet Flow”, .4/4.4 paper 93-0653, NASA TM 105985, 1993. 

Hayder, M. E., Zhou, Y. and Rubinstein, R., “Numerical Simulation of the Mixing Noise in 
Turbulent Flows' , Proceedings of the Fluid Engineering Division Summer Meeting , FED-Vol. 
238, ASME, pp 479 - 484, 1996. 

Hu, F. Q., “On absorbing boundary conditions for linearized Euler equations by a Perfectly 
Matched Layer”, Journal of Computational Physics, 129: 201-219, 1996. 


Hu, F. Q., Hussaini, M. Y. and Manthey, J. L., “Low-Dissipation and Low-Dispersion Runge- 



Kutta Schemes for Computational Acoustics’' Journal of Computational Physics , 124: 177- 
191, 1996. 

Mankbadi. R. R., Hayder, M. E. and Povinelli, L. A., “The Structure of Supersonic Jet Flow 
and Its Radiated Sound" AIAA Journal , 32: 897-906, 1994. 


Ta'asan, S. and Nark, D. M., “An Absorbing Buffer Zone Technique for Acoustic Wave 
P ropagat ion" , -4 1 A A paper 95-0 164, 1 995. 




Figure 2: A snapshot of pressure 
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Figure 8: Free shear layer 







Figure 11: A snapshot of pressure in the jet 



Figure 12: Axisymmetric jet (Linearized Euler) 
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